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Abstract 

Using a powerful combination of Coleman's instanton technique and the method of Banks and 
Bender, the exponential factor for the zero temperature rate of tunneling out of metastable vacuum 
in a system of two identical capacitively coupled phase qubits is calculated in closed form to second 
order in asymmetry parameter for a special case of intermediate coupling C = Cj/2. 
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Introduction 





Tony Leggett's beautiful insights P3EJ into the nature of macroscopic quantum phenomena have long 
been an inspiration to many students of quantum mechanics. One particularly interesting system 
whose macroscopic variables exhibit quantum mechanical behavior predicted by Leggett is the current 
biased Josephson junction (CBJJ). Since the mid 80's a considerable amount of experimental and 
theoretical work has been done to understand its properties It was eventually realized that 
due to its characteristic energy spectrum a CBJJ may play an important role in the field of quantum 
computation 0j. The lowest quasistationary energy levels of the junction, whose spacings are relatively 
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easy to manipulate, may be viewed as the levels of a qubit, and if several such qubits are joined together 
a quantum computer could be built and operated. 

The physics of an isolated junction is by now well understood. However, when two or more such 
junctions are coupled together the calculations of their properties become much more involved. This 
is because we have to deal with higher dimensional problems. The matter is further complicated by 
the presence of additional velocity dependent terms in the corresponding Lagrangians. 

One important property that needs to be addressed is the rate of tunneling out of metastable 
vacuum in a system of several junctions. Such tunneling process (appearance of a voltage drop 
across the junctions) is often used to perform a readout on qubits, and therefore a good theoretical 
understanding of it is essential to the general understanding of CBJJ-based quantum computers. 

In this paper we propose a method for calculating the zero temperature tunneling exponential 
factor. It combines Coleman's instanton technique ,5. and the approach of Banks and Bender [H] 
based on Jacobi's variational procedure. 

With respect to the proposed method the following must be kept in mind: 

1. The tunneling rate is calculated using Coleman's formula (|52|). The action is for the path of zero 
energy which starts at the top of the inverted potential Ue (corresponding to the metastable 
minimum of U), approaches the turning point of Ue (the exit point of U), then goes back to the 
top. Such path is called a bounce. 

2. The action of the system is stationary under variations of the bounce. Being a path of zero 
energy, the bounce is therefore the most probable escape path for tunneling out of metastable 
vacuum. This indicates that the results obtained here should be in agreement with those derived 
using the WKB approximation. 

3. In the presence of velocity dependent couplings, even for highly symmetrical potentials (double 
wells, etc.) the bounce may not be a straight line. In general its shape depends on the interplay 
between the various kinetic energy terms. 

4. Because the action appears in the exponent, the calculated tunneling rate strongly depends on 
how much the bounce actually curves. 

5. In our problem, which deals with identical qubits, only one metastable well is considered. The 
relevant bounce at zero bias difference is indeed a straight segment despite the presence of a 
velocity dependent term in the Lagrangian. 



We are interested in a system of two phase qubits coupled by a capacitor [7j as depicted on Figure 
1. The bias currents are denoted by I\ and I2, C\ 2 are junctions' effective capacitances, C is the 
coupling. Considering the balance of currents in the circuit 



and using Josephson's equations for the voltage and phase difference across a junction whose critical 
current is Iq, 
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Lagrangian for capacitively coupled qubits 



h = I j, + Id + Ic, h = Ij 2 + Ic 2 - Ic 




I j = Iq sin <f>, (f> = aV, 
we are led to the equations of motion for the two superconducting phases 



(2) 












(3) 
(4) 
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where a = h/2e. The corresponding Lagrangian is given by 



L = a 2 



(Ci + C)<j>i (C 2 + C)4> 2 2 ■ • 

o 1 o C0102 



where 



U(^i,4>2;h,h) = -a (% cos0i + Ji0i + J 2 cos 2 + hfa) 



(5) 



(6) 



and the dot ' represents differentiation w.r.t. the time t. 

Denote ji = Ii/Ioi- As ji t 2 — * 1~, 0i, 2 — > n/2. In this limit, 

U ~ - (eiji + e 2 j 2 )7r/2+ ei(l - ji) (^>i -7r/2) + e 2 (l -ja) (^2 
- ei/6 (0i - tt/2) 3 - e 2 /6 (0 2 - vr/2) 3 , 



7T/2) 



(7) 



where e, = a/oi are junctions' Josephson energies. Shifting the local energy minimum to the origin 
gives 



U 



1 r 
6 

27 

T 



ei0 x (a - 0i) + e 2 2 (& ~ 02 ) 
2 



where 



a = 3^/2(1 -ja), 6 = 3^/2(1 -j 2 ) 
are the exit points of the potential in the directions of the coordinate axes and 



(8) 
(9) 

(10) 



are the corresponding potential heights. 

The use of instanton technique [H] for tunneling rate calculations requires the introduction of the 
Euclidean Lagrangian and equations of motion (in non-physical time r), which in this paper will be 
used interchangeably with the usual ones. 

Introducing a small asymmetry parameter rj such that b = (1 +rj)a, the full Euclidean Lagrangian 
becomes 



where 



■ 2 -2 

Le = ~ 1 7, ™30102 - UE(<h,<h\ V) 



Usifaifhw) = -U((t>i,(/)2;v) = \ W<t>\ + £2<t>l ~ o(ei0i + e 2 0f) - r?ae 2 0| 



(11) 

(12) 



rrij = a 2 (Cj + C), and 7713 = a 2 C . 

For identical qubits, C% = C 2 = Cj, /oi = -^02 = -^Cb m\ = ^2 = m + ni3 = a 2 (Cj + C) , m = a 2 Cj, 
e\ = e 2 = e, so that 



Le(4>i,4>2;v) 



m + rri3 



(• 2 ■ 2\ ■ C j 

^1 + 02 J - m 3 0i0 2 - g 0i + 2 - ct(0i + 2 ) - r?a0 2 



We will also need an expression for in the "diagonal" coordinates 

1 1 

xi = -/= (01 + 02) , x 2 = (01 - 2 , 



(13) 



(14) 



which is 



L E (x 1 ,x 2 ;r]) 



mx\ 2 



+ 



Mx? 2 
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2 i 2 

+ x 2 



-rja(x\ 



X 2 



(15) 



where M = m + 2m% = a 2 (Cj + 2C). The three-dimensional graphs of the potential in symmetric 
(77 = 0) and asymmetric {q = 0.30) cases are given on Figures 2 and 3. Here, m = 1 and e = 1. 
In the xi^-coordinates the critical points of Ue are (0, 0) - local maximum, (^/2a/3, \/2 a /3^ and 

((l + r))V2a/3, -(1 + rj)V2a/3\ - saddle points, and ((2 + rj)V2a/3, -r/^a/3^ - local minimum. 
Figures 4 and 5 demonstrate the locations of these points on the corresponding contour diagrams in 
the symmetric and asymmetric cases. 

It is important to clearly understand the meaning of the asymmetry parameter rj. It is due entirely 
to the difference in the bias currents ji . The reference current is arbitrarily chosen to be ji . Experiment 
begins with both bias currents pre-set at some j'imtiai = ji- This corresponds to a given value of the 
exit point a in the (pi direction. Any change in the tunneling rate is always measured in relation to 
the a so chosen. Varying j 2 changes the effective width and the height of the potential barrier in the 
4>2 direction. An increase Aj in j 2 effectively reduces the size of the barrier, resulting in a negative 
value of 7], which is explicitly given by 



Aj 



1. 



1 jinitial 

That, in turn, leads to higher tunneling rate as will be rigorously demonstrated below. 



(16) 



3 A note on instantons 

In the symmetric case {rj = 0), the equations of motion derived from (|13|) and (|15|) using the usual 
variational procedure are 



(m + m 3 )0i - m 3 (f)2 
(m + m 3 )(f)2 ~ m 3 0i 



-~ (30? - 20(h) 
-|(3^-2a^ 2 ) 



(17) 



and 



xi 
x 2 



e 
6m 
e 



2axi 



6M L\/2 
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-X\x 2 — 2ax 2 



(18) 



We now ask: What is the effect of coupling C ~ m.3 on the instanton solutions of ()17|) satisfying the 
boundary and energy (E = 0) conditions 



cp(±oo) = 0, 0(0) = 0? 

When C = (decoupled qubits), there are three such solutions: 
1. one passing through the saddle point of Ue = — U in the <pi direction, 



(19) 



cosh 2 (wr/2)' 



<4V) = 0, 



(20) 



4 



2. another passing through the saddle point in the 02 direction, 



4>f\r)=0, 4V)= ° (21) 

cosh {lot/2) 

and 

3. the third one in the x\ ~ fa + 02 direction of JTe's local minimum (C/'s local maximum), 

4° } (r) = — ^-r, 4V) = 0. (22) 
cosh {lot/2) 

Here, the superscript (°) reminds us that 7/ = 0. Also, 



w - 1 d2U - (23) 

is the frequency of small normal oscillations in the x\ direction. This frequency is also called the 
attempt frequency because it appears in the tunneling rate prefactor. It can be shown by direct 
calculation that for a slightly asymmetric U the attempt frequency acquires a linear in n contribution, 
to — ► lo(1 + 77/4). 

Notice that in configuration space all three instantons are straight segments. Also notice that (|22[) 
should be regarded as non-physical since the two qubits are completely decoupled in this case |S]. 

When C > 0, there are no longer instanton solutions in the form of straight segments passing 
through the saddle points. Indeed, if, for example, we choose to look for a solution in the form 



0, 02 = Mr), (24) 



we would get the system 



"7-302 = 0. 

e 



(m + m 3 )0 2 = -1 (302-2a0 2 ) , (25) 

with two constant solutions 02 = and 02 = 2a/3, which are not instantons. So (|2U|) and (|21|) do 
change when C > 0. 

In this regard we notice that, according to (|17|) and 1)191) . in the extreme case of very strong 
coupling (C » Cj), 0i (t) - 02 (t) ps =>- 0i(t) - 2 (r) w At + -B ^ x (r) w 2 (r), so both 
instantons disappear. Apparently the velocity dependent coupling becomes so strong that it prevents 
the formation of the corresponding bounces. 

However the xi-instanton described by 1)22(1 becomes physical, surviving any coupling. Any change 
in it is due to asymmetry in the potential. 

We now come to a subtle point: Do the 0i- and fa-bounces actually exist at finite coupling C > 0? 

Here is an argument I would like to put forward regarding this matter: Assume all three exist. 
Then the tunneling rate is roughly the sum of three contributions, e -5 ^!^ + e -1 ^/^ + e~ Sx t h , where 
various S"s are the actions for the corresponding instantons. Now, as C — > 00, P — > e~ Sx ^^ h , which is 
equivalent to 

fW^oo. (26) 



>4>i ' 4>2 

However, as we have just shown, in this limit the bounces must satisfy the condition 0i(r) ~ 02 (t). 
This contradicts (|2T)|) because in the fa ps 02 > region U is continuous and bounded. We therefore 
conclude that at any finite coupling there is one and only one instanton given by (|22|). 
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The physical meaning of this conclusion is the following: When two qubits are capacitively coupled 
there can be no "tunneling" voltage drop on one of them without having a comparable drop on the 
other. Apparently the qubits always "tunnel" together. This is especially obvious in the limit C —* oo 
when the qubits effectively become one. 



4 Jacobi's variational principle 

We now turn to the effect of asymmetry (rj ^ 0) on the instanton solution (|22j). It is convenient to 
work with the equations of motion (with energy condition E = 0) derived via Jacobi's variational 
principle 

5 j ds^2U{s) = (27) 

with the constraint 

ds 2 = m(dxi) 2 + M(dx 2 ) 2 . (28) 
Performing the variation we get (c/. [0]) 

d 2 x\ dx\ ( dx\ dU dx 2 dU \ dU 
2Um - n + m— — — — — h 



ds 2 ds V ds dx\ ds dx 2 ) dx\ ' 

2UM^ + M^(^^ + ^%L) = *L. (29) 
as z ds V ds ox i ds 0x2/ 0x2 

Let us assume that these equations have been solved in the symmetric case (77 = 0): x(°)(s) = 
(x^\s), x^2°\s)). The solution in the asymmetric case (77 7^ 0) will be written in the form of a power 
series in rj: 

x(a) = x(°) (s) + jyxW (s) + if^ (s) + . . . . (30) 

From (j28j) several important relations among various x.( k \s) can be derived. [From here on, the 
dot ' stands for the derivative with respect to s and the dot product is evaluated with the help of the 
metric diag(m, M).\ To second order in rj we have: 



1 = x • x = (x(°) + ryxW + r/ 2 x( 2 )) • (x(°) + rpcW + 

= l + 27 ? x( ).x( 1 )+r ? 2 (2x( )-x( 2 )+xW-x( 1 ) N ). (31) 



This leads to 

xW • i« = (32) 

and 

x(°)-x( 2 ) = -ix( 1 ).x( 1 ). (33) 

We will take x^ )(s) to be the s-parametrized form of the instanton solution x(°)(t) given by (|22|) . 
Thus in configuration space x^(s) is a straight segment 

x (0) (s)=s(l/Vm, 0), 0<s<si, si = aV2m. (34) 

Now, from (|32l) and (|34|) it follows that the first component of x^ 1 ) (s) must be zero which allows 
us to choose x^ 1 ) (s) in the form 

x«( S ) = (0, X (s)). (35) 
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Thus to second order in rj, 

x(s) = (xi(s), x 2 (s)) = (xf\s) + ri 2 xf\s), 7]\{s) +rj 1 xf ) (t 
and the potential is given by 

-M(r,)/e = (xS°)) 2 (i=4 0) -a) 

-J7|(4 0) ) +V 2 "^ x i 0) (4°M 2) +X 2 ) + a (x^ 0) x - 2x[ 0) a;f ) - x' 



(36) 



mv/2m 



(s - si) 



2 2 
7/S Si 7] 

m\/2m 2 y^^m 
Prom (J231), (EH), and (jSHJ) it follows that 



3s 



S X 



(2) 



± + X 2 )+s 1 



s ( x ~ 2xf } 



m 



1 M 



rn 



X 



(37) 



or 



4 2) (s) 



2^ X ' 
1 M 



2 v/m 



ds x (s) . 



(38) 
(39) 



We see that, as pointed out in jS], the calculation of tunneling effects to second order in r\ requires the 
knowledge of the tunneling path to only first order in rj. 
Thus, to first order, 



dU(s) 


es 


dx\ 




dU(s) 


erj 


dx 2 


6y/2m 



(3s - 2si - r]si) , 
2(3s-si) X + ^) , 



and from l)29j) we get the following equation for f(s) = y / mx(s): 

2s 2 (l - s/s x )/(s) + s(2 - 3s/ Sl )/(s) - 2k(1 - 3s/ Sl )/(s) 

where 

K = m/M = Cj/(Cj + 2C) 

characterizes the strength of the capacitive coupling. 
We may further perform the change of variables 

1 /2 

/0) -> g{z) = f(z)/si, s -> z = (1 - s/si) 1 , 
with z ranging from 1 to 0, transforming (|42|) into 



-ks, 



(i-^) 



2^ A^) 2z- d5 '^' ) 



dz 2 



dz 



+ K 



12 



1-z 2 



9{?) 



-2k. 



(40) 
(41) 

(42) 
(43) 

(44) 
(45) 



This we recognize as an associated inhomogeneous Legendre equation whose homogeneous solutions 
are associated Legendre functions Pjf and Q 1 ^ of degree v (where v{y + 1) = 12k) and order fj, = 2yfR. 
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5 The case of intermediate coupling C = Cj/2 

Up to this point the derivation was valid for any coupling. However equation (|45jl is difficult to work 
with in general. To simplify the presentation we consider a special case k = 1/2 which corresponds to 
intermediate coupling C = Cj/2, 



(1-*) 
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l-z 2 



dz 2 dz 

The solution to (|46j) satisfying the all-important initial condition g(l) = is 

g{z) = {z 2 - l)/2 

which gives 



1. 



X{s) 



x(o) = o. 



This leads to a very simple expression for the potential U in (j37j) . 

,2 



6m\/2m 



(^exit s) ) 



where 



is its exit point. 



Scxit = si ( 1 + -?? + 



(46) 

(47) 
(48) 

(49) 
(50) 



6 Results of numerical simulations 

A Runge-Kutta simulation of the xi-instanton trajectory in configuration space is presented in the 
following table (where, again, r] = 0.30, m = 1, and e = 1): 



s 


„sim. 
X 1 


prcd. 

x 1 


„sim. 

x 2 


pred. 
x 2 


0.000 


0.000 


0.000 


-0.000 


-0.000 


0.500 


0.491 


0.489 


-0.066 


-0.075 


1.000 


0.983 


0.978 


-0.130 


-0.150 


1.500 


1.475 


1.466 


-0.193 


-0.225 


1.654 


1.626 


1.617 


-0.213 


-0.248 



(51) 



The predicted values of x^ 6 ' and x^ e ' are given to second and first order in n, respectively, since 
that is the accuracy to which they are needed for the second order tunneling rate calculations. 

The corresponding exit point has also been found numerically. It is s|^' = 1.654, which is in good 
agreement with the value s^i? = 1-658 predicted using (|5U|) . 



7 The exponential factor 

The tunneling exponential factor P can now be calculated via 

P = exp (-S Q /ft), (52) 
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where 

S = 2 J ds yJ2U(s) (53) 

is the action for the tunneling path x(s). The factor of 2 is due to the bounce. 
Using we find: 




where the terms to second order have been retained and r\ is given by (|16|). 

Notice that the zeroth order term (corresponding to n = 0, rj = 0; strongly coupled qubits, zero 
bias difference) coincides with the well known result of Caldeira and Leggett [2] (for instanton of 
amplitude ay/2). Also, negative r\ (higher bias) increases the tunneling rate, as expected. 

8 Conclusions and the outlook 

In this paper a combination of instanton technique and Jacobi's variational procedure has been used for 
calculation of the rate of tunneling out of metastable vacuum in a system consisting of two capacitively 
coupled phase qubits. In the interesting case of intermediate coupling the tunneling rate was explicitly 
calculated in closed form to second order in asymmetry parameter rj. 

As a by-product, a detailed analysis of instanton trajectories in the presence of velocity dependent 
coupling has been performed. It was found that a curved version of only one such trajectory, (|22jl. 
gives the major contribution to the tunneling rate in a capacitively coupled system. Physically this 
means that simultaneous tunneling of both qubits is overwhelmingly more probable than any other 
mode of tunneling. 

Tunneling out of higher energy levels and tunneling out of various superpositions of the ground 
and excited states is an active area of modern experimental research. While the methods of this paper 
can be easily applied to tunneling out of metastable vacuum, it is important to know whether more 
complicated tunneling processes can be analysed using the same approach. 

The author believes that generalization to higher energy states is straightforward. Again, the most 
probable escape path would be given by a curved version of the bounce along the x\ axis, similar to 
1|22|). This time, however, the bounce would pass through the exit points of Ue corresponding to non- 
zero energy. The rate of tunneling out of a superposition then would be regarded as an uncorrelated 
sum of tunneling rates out of sharp energy eigenstates multiplied by the corresponding weighting 
factors. 
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